Reentrant superconductivity in a strong applied field within the tight— binding model 
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It was suggested by Rasolt and Tesanovic that the Landau level quantization in isotropic super- 
conductors could enhance superconductivity in a very strong magnetic field, above the upper critical 
field (H C 2). We derive a generalized Harper equation for superconducting systems, and show that 
a similar reentrant behavior appears in a lattice model, even though the Landau-level-structure is 
destroyed by the periodic potential in that case. Both the orbital and the Zeeman field-induced 
effects are taken into account. 
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I. INTRODUCTION 

There are two mechanisms responsible for a suppres- 
sion of conventional superconductivity in an external 
magnetic fieldEl: the Pauli pair breaking and the dia- 
magnetic pair breaking. First of them, the Pauli pair 
breaking, is connected with the Zeeman coupling. The 
magnetic field tends to align the spins of the electrons 
forming the Cooper pair, and the singlet superconduc- 
tivity, disappears at the Chandrasekhar-Clogston (CC) 
limitS. However, this critical field for the majority of 
type-II systems is found to be above H C 2 determined by 
the orbital (diamagnetic) pair-breaking. Especially, this 
effect is of minor significance in materials with low effec- 
tive g-factor. Another possibility is the superconductiv- 
ity with nonhomogeneous order.-parameter (the Larkin- 
Ovchinnikov-Fulde-Ferrell statea), which can exist above 
the CC limit. One can also look for high-magnetic-field 
superconductivity in superconductors with triplet equal 
spin pairing. 

The second effect, the diamagnetic pair-breaking, usu- 
ally crucial in determining the upper critical field, is con- 
nected with the orbital frustration of the superconduct- 
ing order parameter in a magnetic field. This frustration 
enlarges the free energy of the superconducting state, 
and, when the magnetic field is strong enough, the nor- 
mal state becomes energetically favorable. The orbital 
effect can be reduced in layered two-dimensional super- 
conductors, when the applied magnetic field is parallel 
to the conducting iayers. Such a situation has been ana- 
lyzed theoreticallyu and recently observed experimentally 
in organic conductorsu. 

However, it was shown that large values of the 
critical field are possible also in the systems without 
two-dimensional layers, i.e., in systems where the or- 
bital effects are present. When describing a supercon- 
ductor,-|Within the Ginzburg-Landau-Abrikosov-Gor'kov 
theoryQ, one treats the magnetic field in the semiclassical 
phase-integral approximation, thus neglecting the quan- 
tum effects of the magnetic field. This approximation is 
valid for relatively small fields, when huj c <C ksT c (or 
lu c <C 2tt/t for large impurity concentration, where r is 
the elastic scattering time). In this regime, the number 
of occupied Landau levels is very large and the energy 
spacing of them is very small, and therefore this discrete 



structure is not observable. However, when the mag- 
netic field increases, the Landau-level degeneracy also in- 
creases, so the number of occupied levels decreases and 
one has to take it into account. The inclusion of the Lan- 
dau level quantization in the BCS theory leads to reen- 
trant behavior at a very high magnetic field (huj c 3> ep)- 
Namely, when only the lowest Landau level is occupied, 
T c is increasing function of H, limited only by impurity 
scattering and the Pauli pair breaking effect. 

The aim of this paper is to show that the reentrant 
behavior survives in the presence of a strong periodic 
lattice potential. 

Weak, unidirectional periodic potential removes (or, at 
least, modify) the Landau-level structure: the levels are 
broadened (they form "Landau bands") and the degen- 
eracy is lifteduJ. The width of a Landau band oscillates 
as the magnetic field is tuned as a consequence of com- 
mensurability between the cyclotron diameter and the 
period of the potential. It results in a magnetoresistance 
oscillations (Weiss oscillationsO) . If the periodic poten- 
tial is modulated in two dimensions, "minigaps" open in 
the "Landau bands" , and the energy spectrum of the sys- 
tem plotted versus the. applied field composes the famous 
Hofstadter butterflyuEl, recently observed experimentally 
in the quantized-Hall-conductance measurementE3. The 
same spectrum can be obtained in a complementary 
limit, when the lattice potential is strong (tight bind- 
ing approach) and the field is weak. It is interesting that 
when the periodic potential does not lead to a scattering 
between states from different Landau levels, the eigen- 
value equations in both the limiting cases are formally 
the sameo. Of course the parameters have a different 
physical meaning. 

The simplest model for the case where a aplied field 
and a lattice potential are present simultaneously, is com- 
monly|-|Ceferred to as the Hofstadter or Azbel-Hofstadtcr 
modelBU. The corresponding Hamiltonian describes 
electrons on a two-dimensional square lattice with the 
nearest-neighbor hopping, in a perpendicular uniform 
magnetic field. The Schrodinger equation takes the form 
of a one-dimensional difference equation, known _as_jtbe 
Harper equation (or the almost Mathieu equation )l§E3Li± 
It is also a model for one-dimensional electronic sys- 
tem in two incommensurate periodic potentials. The 
Harper equation has also links to many other areas of 
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interest, e.g., the quantum Hall effeC|L-|quasicrystals, 
localization-delocalization phenomenaliJ't£l, the .eoiicom- 
mutative geometry^, the renormalization groupOlla, the 
theory oiJractals, the number theory, and the functional 
analysiscj. 

The Hofstadter model is useful in approach to the fun- 
damental problem of the external magnetic field influ- 
ence on the superconductivity. Most of the works de- 
voted to superconductors in the mixed gtjate are based 
on the Bogolubov-de Gennes equationso, particularly 
useful for spatially inhomogeneous syjstems, e.g., for an 
isolated vortexcll or a vortex latticeE3. However, in the 
regime H c i <C H <C H c2 we can neglect contributions 
to the spectrum from the inside of the vortex core (for 
H <C H C 2 the distance between the vortices is large) and 
regard the magnetic field as uniform in the sample (for 
H c i <C H). We derive, under these assumptions, a lattice 
model for the superconductor in applied field (in the nor- 
mal state such a system is described by the Hofstadter 
model). In this paper we present a generalized Harper 
equation that describes the influence of magnetic field 
on the two-dimensional tightly bound electrons in the 
superconducting state. 



II. THE MODEL 

In analogy to Hofstadter's approach, we couple 
the magnetic field to the system via the Peierls 
substitutions, i.e., multiply the hopping matrix elements 
by a phase factor which depends on the field and on the 
position within the lattice. Thus, the vector-potential- 
dependent hopping integral for sites i and j is given by 



hj (A) = iexp ( — 



R, 



A-tfl 



where t is the usual hopping integral. We also include 
the Zeeman term. In effect, the BCS Hamiltonian has 
the form 

(ij),<j i,tr 



(2) 



(ij) 



where the Zeeman splitting is given by e a = — ^ag^sH, 
(7 = 1 for spin up and —1 for spin down, g is the 
Lande factor, and [i is the chemical potential. Here, we 
have introduced the spin-singlet pair amplitude Ay = 
\(ci\Cji — CiiCji). The strength of the nearest neighbor 
attraction V is assumed to be field independent. The va- 
lidity of this assumption depends on the nature p£ pairing 
potential and the strength of the magnetic fieldEa For ex- 
ample, in the t — J model Jy(A) = 4tij(A)tji(A)/U is 
strictly field independent, since the change of the phase 
generated when an electron hops from site i to j and 
back, cancels out. Such an assumption has been also 



partially justified on the basis of antifecromagnetic-spin- 
fluctuation-driven superconductivity.t3 

Our starting point is a two-dimensional square lattice 
with basis vectors a = (a, 0,0) and b = (0, a, 0), im- 
mersed in a perpendicular, uniform magnetic field H = 
(0,0, H). We choose the Landau gauge, A = (0, Hx,0). 
Since the vector potential A is linear in x, the transla- 
tion corresponding to the vector a shifts the phase of the 
wave function. This shift can be compensated by a gauge 
transformation, introducing magnetic translations. If the 
magnetic flux per unit cell, is a rational multiple of 
the flux quantum $o = he/ e, i.e., if 
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with p and q coprime integers, we can define magnetic 
lattice, with q& and b as the basis of the magnetic unit 
cell. Such an enlarged unit cell is penetrated by p flux 
quanta. Magnetic translations corresponding to the mag- 
netic lattice vectors (R = nqa+mh, with n, m - integers) 
commute with each other and with the Hamiltonian. If 
the system is of rectangular shape with L x sites in the x 
direction and L v sites in the y direction, and L x is a mul- 
tiple of q, we can find eigenfunctions which diagonalize 
the Hamiltonian and the magnetic translation operators 
simultaneously. Due to the absence of translational in- 
variance with vectors mb, vectors k = (k x , k y ) from the 
firs Brillouin zone (\k x \ < n/a, \k y \ < 7r/a) are not good 
quantum numbers. Instead, we have to use vectors from 
a magnetic (reduced) Brillouin zone (MBZ), defined by 
\k x \ < Tt/qa, \k y \ < n/a, to enumerate the eigenstates. 
The Hamiltonian (1) in the momentum space can be writ- 
ten as 



(1) H = -ty] b cos (k x a) c£. -CT c k , CT 



k,<r 

, -ikyu f , ikya t „ 

°k-g,cr t ' k ! CT t e L k+g,<T Lk <°' 

kcr k 



H.C 



■()•) 



where 



k' 



Vk,k'(c_k'l c k'T)7 



(5) 



and g = (2irp/q, 0). In order to rewrite the above Hamil- 
tonian as a sum over the MBZ we introduce a multicom- 
ponent Nambu spinors 



~ ( c k, T ' c k-g,T> c k-2 g . j- 



C~k,l7 c -k+g,l; C_k+2g,J.j 

Then Eq. (4) can be written as 
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H — ^2 Cl_H k C k , 
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where the prime denotes summation over the MBZ and 
has a block structure 







A* 


-T-ki 



(8) 



The diagonal blocks describe noninteracting lattice 
fermions under the influence of magnetic field, andphave 
the form similar to that derived by Hasegawa et alEB 
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(9) 

where M„ j(T = 2 cos (fc x a + n-f) + e a - fi, 7 = |g| = 2np/q, 
and diagonal matrix Ak represents the pairing ampli- 
tudes 



A k = diag(A k , A k - 



g' 



A 



k-( 9 -l)gJ 



(10) 



Diagonalization of Eq. (8) provides a set of eigenenergies 
{£k,i}, where i enumerates 2q values corresponding to a 
given k from the MBZ. 

The pairing amplitude in the presence of external mag- 
netic field is determined self-consistently from the BCS- 
like equation 
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where N — L x L y and the prime summation denotes 
again summation over the MBZ. In the following we re- 
strict ourselves to the singlet pairing in the s-wave chan- 
nel (Ak = A), even though Eq. (11) is completely gen- 
eral. Generally, this equation can be used, e.g., to ana- 
lyze the magnetic-field-induced change of gap parameter 
symmetrycJ or the upper critical field in the systems with 
the spin-triplet pairing. The latter area of application is 
especially attractive, since in these systems the Pauli pair 
breaking mechanism is absent, and the upper critical field 
is expected to be very high.EEl 



III. RESULTS 
A. Orbital effects 

The transition lines T C (H), obtained in the absence of 
the Zeeman splitting (g = 0), are presented in Fig. 1. 

Note that for a of the order of a few angstroms, experi- 
mentally available magnetic fluxes are much less than <& . 
Consequently, these plots correspond to the region of ex- 
tremely high magnetic field. The size of the Hamiltonian 
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FIG. 1: Critical temperature T c as a function of $/$o. Upper, 
middle, and lower panel shows results for V/t = 7, 2 and 1, 
respectively. Small arrows indicate the transition temperature 
calculated from the usual BCS equations in the absence of 
magnetic field. 
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FIG. 2: Critical temperature T c as a function of $/<hriin the 
low field regime for V/t =1. Results taken from Ref.E3 



matrix H^, that has to be diagonized for all values of k in 
each step of the iterative procedure, is 2q x 2q. Therefore, 
since the magnetic flux is proportional to q^ 1 , the pro- 
posed approach does not allow to carry out calculations 
for a small magnetic field. This is why the transition 
lines in Fig. 1. start at finite magnetic field. 

For weak field thermal smearing and/or disorder in- 
duced broadening destroy the Hofstadter butterfly struc- 
ture. In the absence of the lattice periodic potential this 
regime corresponds to a classical limit, where the num- 
ber of occupied Landau levels is huge, and the Ginzburg- 
Landau description of the mixed state is valid. In this 
regime, in accordance with the common feeling, super- 
conductivity in the tight-binding systpm_is suppressed 
by magnetic field, disappearing at i^Eijtij. The corre- 
sponding transition line is presented in Fig. 2. 

The method used in Ref.Ba does not work at low tem- 
perature and the present method does not work at weak 
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FIG. 3: Normal state density of states for different values of 
magnetic field. The dotted line in the first panel represents 
density of states in a superconducting state for different ratio 
p/q and for A = 0.7t. 
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FIG. 4: Normal state density of states for p/q = 1/2 (a) and 
p/q — 10/21 (b). (c) Comparison of integrated densities of 
states for p/q = 1/2 (dotted line) and p/q = 10/21 (solid 
line) . 



field. Therefore, there is no crossover line from the low 
to high field regimes. 

The transition lines for 1/2 < $/$() < 1 can be ob- 
tained reflecting the lines presented in Fig. 1. around the 
line $/$o = 1/2, and T C (<I>) is periodic on <5> . Both these 
properties reflect properties of the Hofstadter butterfly. 
Of course, these unphysical results are valid only when 
the Pauli pair breaking is neglected. The influence of the 
Zeeman splitting will be discussed later. For strong pair- 
ing potential, comparable with the band width, the crit- 
ical temperature in the reentrant regime is almost field 
independent (see Fig la.). As V is reduced, the influence 
of the nontrivial density r-pf states becomes apparent. It 
was shown by Hofstadtera, that in normal state the Bloch 
band for &/Qq = p/q is symmetric and broken up into q 
distinct energy bands. In the half-filling case the Fermi 
level (Ep) is located in the center of the (unperturbed) 
subband. Therefore, if q is odd, Fp, points to the singu- 
larity of the central subband (a remnant of the original 
van Hove singularity), whereas for even q it is in the gap 
between two subbands (In fact, for even q these subbands 
touch itself at the Fermi level). This is depicted in Fig. 
3. 

The changes of the density of states result in an os- 
cillatory behavior of T C (H): T c approaches its maxima 
for odd q and is reduced for even q. Similar oscillations 
were predicted by Rasolt and Tesanovictl in a homoge- 
neous system, where the Hofstadter spectrum is replaced 
by the Landau-level ladder. 

The superconductivity suppression is especially appar- 
ent for small and even q, when V is comparable with the 
central-gap width. The smooth character of the function 
T c ($) close to ^/^o = p/q and for small q (e.g., close to 
the values p/q = 1/2, 1/3, 1/4), seems counterintuitive, 
since a tiny detuning of the magnetic field completely 
changes the spectrum. For p/q = 1/2 the spectrum con- 
sists of two subbands, whereas for p/q — 10/21 there are 
21 narrow subbands (see Fig. 4). However, in spite of this 
difference, the integrated densities of states, presented in 
Fig. 4c, are almost the same. 

For larger q, the differences between 1/q and l/(q + 
1) are smaller, and consequently the distances between 
successive minima in the density of states decrease. For 



a strong pairing potential (and high T c ) there is large 
number of subbands within a range of energy ~ ^T c 
and then the amplitude of oscillations is strongly reduced. 
On the other hand, for the weak potential (i.e., at low 
temperature), this irregular oscillations are visible even 
at low fields (cf. Fig. lc). 



B. Zeeman splitting 

The previous discussion ignored the effect of Pauli pair 
breaking. We consider it next. Since the Zeeman split- 
ting is proportional to the magnitude of the magnetic 
field and the orbital effect depends on the flux, we have 
to find a relation between these two quantities. It can 
be done by using the relation t = h 2 /2m*a 2 , where m* 
is the effective mass. Then the ZeemaajSplitting is given 
by g f i B H = 2ng*Zt, where g* = g^H 

The inclusion of the Zeeman term results in a reduc- 
tion of the phase space available for pairing. For strong 
pairing potential, when the structure of the Hofstadter 
butterfly is hidden, it leads to a monotonic reduction of 
T c with increasing magnetic field. Such a situation is 
presented in Fig. 5a. 

However, for smaller values of V, when fcgT c is com- 
parable with the miniband (or minigap) widths, the sit- 
uation is more complicated. The Zeeman term leads 
to the splitting of each of the minibands into a spin- 
up and spin-down minibands. To have nonzero T c we 
need minibands of both types present close to the Fermi 
level. As the magnitude of the splitting is proportional 
to the magnetic field, T c will be an oscillatory function 
of the magnetic field. When the spin-up and spin-down 
minibands overlap at the Fermi level, T c is strongly en- 
hanced. This mechanism may induce superconductivity 
in regions, where T c is zero or close to zero in the absence 
of the Zeeman splitting (compare the solid and dashed 
lines in Fig. 5b). For example, for < f>/ < E ) o = 1/4 Ep is 
located in the central minigap for g* — 0, whereas there 
is a singularity at Ep for g* = 0.15. The corresponding 
densities of states are presented in Fig. 6. 
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FIG. 5: Critical temperature T c as a function of <E>/$o in the 
presence of the Zeeman splitting for V/t — 7 (a) and V/t = 2 
(b). 




FIG. 6: Normal state densities of states for $/$o = 1/4. 
Solid and dashed lines correspond to g* — and g* — 0.15, 
respectively. 



IV. DISCUSSION 

Let us comment on the possibility of observing the os- 
cillatory behavior of T c in real systems. Assuming lat- 
tice constant a= 2A the magnetic field required to ob- 
tain <f>/$o ~ 1 is O(10 5 )T, which is obviously too large. 
However, there are some possibilities to overcame this 
problem. For example, it was recently showrtll, that 
in some three-dimensional systems fractal spectra, like 
Hofstadter's butterfly, can be obtained for $/<£>o -C 1. 
On the other hand, it is possible to reach the needed 
increase of flux enlarging the lattice constant. Two- 
dimensional superconducting wire network can be suit- 
able for this task, since the magnetic field corresponding 
to $o is about 1 mT for a network cell of 1 /jm 2 , and the 
system can be mapped onto a tight-binding one. An- 
other possibility is connected with the case, where the 
influence of the modulation potential on the Landau- 
quantized 2D electron system may be considered as a 
small perturbation. This situation is complementary to 
the tight-binding case, but the energy spectrum is also 
obtained by solving the Harper equation. Therefore, one 
can expect similar behavior of T c in 2D superconduct- 
ing systems modulated in two dimensions. Again, since 
the modulation period is larger than the lattice constant, 
the required values of are well within the experimental 
accessibility. 
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